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Abstract. Consider a broken geodesies a([0, 1]) on a compact Rie- 
mannian manifold (M,g) with boundary of dimension n > 3. The 
broken geodesies are unions of two geodesies with the property that 
they have a common end point. Assume that for every broken geo- 
desic a([0,Z]) starting at and ending to the boundary dM we know 
the starting point and direction (a(0), a'(O)), the end point and direc- 
tion (a (I), «'(/)), and the length I. We show that this data determines 
uniquely, up to an isometry, the manifold (M,g). This result has ap- 
plications in inverse problems on very heterogeneous media for situa- 
tions where there are many scattering points in the medium, and arises 
in several applications including geophysics and medical imaging. As 
an example we consider the inverse problem for the radiative transfer 
equation (or the linear transport equation) with a non-constant wave 
speed. Assuming that the scattering kernel is everywhere positive, we 
show that the boundary measurements determine the wave speed inside 
the domain up to an isometry. 

AMS classification: 35J25, 58J45. 
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1. Introduction. 

1.1. Main result. Let us consider a compact Riemannian manifold 
(M,g) with boundary of dimension n > 3. Let SM denote its unit 
tangent bundle. The classical boundary rigidity problem is the fol- 
lowing (see [laillllllllETlEllHaSllETl^]): Assume that we 
know the distances d(x, y) of boundary points x, y G dM. Can we 
determine the isometry type of the manifold (M,g)7 Michel [30l 131] 
observed that in the case of simple manifolds these distance functions 
also determine the values of the bicharacteristic flow at boundary, the 
so-called scattering relation or lens relation, that is, 

C={(x,£), (y, C), t) E SM x SM x E : x,y G dM, 

(7*,e(*), d tlx^if)) = (Vi C) for some t > 0} 

where 7^ is the geodesic of (M, g) that leaves from x to direction £ 
at t = 0. In other words, C gives the information when and where 
and in which direction a geodesic, sent from the boundary, hits again 
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the boundary. It was shown in [16] under some conditions (see also 
[21 Ej) that the wave front set of the scattering operator associated 
to the wave equation for the Laplace-Beltrami operator of a smooth 
Riemannian metric determines the scattering relation. The natural 
conjecture is that for non-trapping manifolds the scattering relation 
determines the isometry type of the manifold. For recent progress on 
this problem see the survey papers [35l 140]. 

In the case of a very heterogeneous media with many scattering 
points inside the manifold one can obtain further information by look- 
ing at the propagation of singularities of waves going through the man- 
ifold. This is the broken scattering relation or broken lens relation that 
we proceed to define. 

A broken geodesic (or, a once broken geodesic) is a path a = a X) ^ z>ri (t), 
where z = J x ,t(s) £ M for some s > 0, rj G S Z M, and 



(See Fig. 1.) In Riemannian geometry broken geodesies are considered 
e.g. in the classical Ambrose theorem [4], which says that the parallel 
translations of the curvature tensor along broken geodesies determine 
uniquely a simply connected Riemannian manifold. 

We denote by £(a x> ^ ZjV ) G M + U {oo} smallest I > such that 
a x£,z,r}(l) £ dM. Denote by v the interior unit normal vector and 



the incoming and outgoing boundary directions respectively. 

The boundary entering and exiting points of broken geodesies define 
the broken scattering relation, 

R = {(x,0,{y,0,t)eSMxSMxR + :{x,OeQ+, (y,0eQ-, 
t = £(a x ^ v ), and 

(ax,z,z,Ti(t),dtoix,z,z,v(t)) = {y,0 for some ( z ,v) e s^"}- 

Our main result is: 

Theorem 1.1. Let (M,g) be a compact Riemannian manifold with 
a non-empty boundary of dimension n > 3. Then dM and the bro- 
ken scattering relation R determines the isometry type of the manifold 
(M, g) uniquely. 

We remark that this result doesn't assume any a-priori condition on 
the metric g or the manifold M. The difficulty in proving the result 
lies in the possible complicated nature of the broken geodesic flow. 
The proof of the theorem above and the other results stated in the 
introduction are given in sections 2-3. 




by 



Q+ = {(x,£) G SM : x G dM, (f , u) g > 0} 
n_ = {(x, £) G SM : x G dM, (f , u) g < 0} 
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Figure 1. Left: Propagation of singularities and mul- 
tiple scattering for the radiative transfer equation. 
Right: A broken geodesic corresponding the relation 

(Oo, fo), (a?i, 60, *) e i? with t = s 1 + s 2 . 



1.2. Application: Radiative transfer equation. As mentioned ear- 
lier the broken scattering relation can be determined by probing with 
waves a very heterogeneous medium with many scattering points and 
observing at the boundary the effects. The strongest singularities of the 
waves are the ones propagating through the medium without any re- 
flection and this determines the scattering relation. The next stronger 
singularities correspond to the waves reflecting only once and this de- 
termines the broken scattering relation at the boundary. This type of 
situation arises in geophysics due to the many discontinuities in the 
surface of the earth that act as reflectors and in optical tomography, 
a novel medical imaging technique that allows one to reconstruct the 
spatial distribution of optical properties of tissues by probing them by 
near-infra-red photons [6l [3, [£3, CEH EU] . This can be formulated as an 
inverse problem for the radiative transfer equation and we consider this 
application in more detail below. For previous mathematical analysis 
on the problem, see e.g. [8l [TOl [TTT |2TI [22l HO] . 

To avoid artificial difficulties on how to formulate the boundary value 
problem for the radiative transfer equation, we consider a non-compact 
complete manifold (A, g) without boundary. The inverse problem we 
study is to find the metric in a compact subset M with smooth bound- 
ary using external measurements made in the set U = A \ M. 

We say that the function u(t, x,£) defined on (t, x, £) G [0, oo) x SN, 
is a solution of the radiative transfer equation on A if 

(1) (Hu)(t, x, + a(x, Ou(t, x, - (Su)(t, x, £) = 0, 

u(t,x,£,)\ t=0 = w(x,£). 

Here H is the bicharacteristic flow on the tangent bundle TA, 
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where (x 1 , . . . , x n , £1, . . . , £ n ) denotes local coordinates on the tangent 
bundle TN corresponding to local coordinates (x 1 , . . . ,x n ) of M and 
= gi k £ k . The operator S, called the scattering operator, is 



Here K G C°°(SN <S> SN) is called the scattering kernel and c n = 
vo\(S n ~~ 1 ). Finally, the function a G C°°(SN) is called the attenuation 
function. We denote the solution of (pQ) with the initial value w G 



For the results concerning the radiative transfer equation we need a 
few more definitions. We say that the complete manifold N is simple if 
for any x,y G iV there is only one geodesic connecting these points. We 
say that M C iV is strictly convex if all points in M can be connected 
with a geodesic segment lying in M and the second fundamental form 
of DM is positive. 

We say that scattering kernel K is positive in M mt if 

K(x, f , O > 0, for all x G M int and f , f G S X N. 

Next we define the external measurements. We assume that for any 
w G C^°(SN), such that w(x, £) = for i 6 M we know solution 
■u™(x, £, t) for x & U. In other words, we assume that we are given the 
measurement map A : C^(SU) -> x 5E/), 



Note that the map A gives us the geodesic flow in U and thus it deter- 
mines the metric gij(x) for x G [/. Also, it can be used to determine 
the absorption a\u- 

Theorem 1.2. Lei N be a complete simple manifold, M C N a 
compact and strictly convex set with smooth boundary. Assume that 
K(x, 9, 9') vanish for x^M, that is, K G (7 °°(SM ® SM) and that K 
is positive in M mt . 

Moreover, assume that we are given the setU — N\M and the mea- 
surement map A. These data determine uniquely the broken scattering 
relation of the manifold (M,g). 



2.1. Auxiliary Lemmata. Let (M,g) be a compact manifold with 
boundary, DM. In the following, we use an auxiliary smooth closed 
compact n-manifold (M,g) that contains (M,g). We continue to use 
notation ■y x ^(t), (x,£) G SM, for the geodesies on M with 7x,§(0) = x 
and 7^,c(£) = £■ All geodesies are parameterized by the arclength. We 

denote by dist^(x, y) and dist(x,y) the distance functions on M and 




C°°(SN) by u(t,x,0 = u w (t,x,0- 



Aw = u w \ R+xS u- 



2. Proof of Theorem 11.11 
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M, respectively. To simplify notations, we denote 

(x ,^ )R t (x 1 ,^ 1 ) ifandonlyif ^(x , £ ), (x u -£i),tj G R. 

On M and M, we will use various critical distances along geodesies. 
We start with critical distances associated with the Riemann exponen- 
tial map, exp^,, 

exp x : T X M = S X M x R + — > M, expM) = T*,^), 

£ G S X M, s G R+. The cut locus distance along 7^, denoted by 
tr(x,£), is defined by 

(2) t r (x,£) = max{s > : dist^(x, 7^(5)) = s}. 

The cut locus distance tr(x, £), (x, £) G «SM determines the injectivity 
radius inj (M) of M, 

inj(M)= mm_T R (x,£). 

(x,t)eSM 

We say that the set 

uj x = {yEM: y = i x &{t r {x, 0), £ e S X M}, 

is the cm£ loc«5 with respect to x. The cut locus ui x consists of two 
types of points. We say that a point y G u x is an ordinary cut locus 
point if there are £, 77 G S^M, 77 7^ £ with 

Tr(z, = r fl (a;, 77), 7^(r fi (a;, 0) = 7x,i,(tr(z, v)) = V- 

Consider now the differential of exp^ at s£ that is denoted by ciexp^ \ s g. 
We say that a point y = 7 x ,g(s) is a conjugate point along 7^, if the 
differential ciexp^, \ S £ : T X M — > T^M is degenerate. This is equivalent 
to the existence of a non-trivial Jacobi field Y(t) along 7 = 7 X ,^([0, s]) 
with the Dirichlet boundary conditions Y(0) = and Y(s) = 0. For 
{x,£) G SM we define the conjugate distance t c (x, £) G M + U {00} to 
be 

t c (x, £) = inf {s > : rfexp^, | S £ is not one-to-one}. 

Each point y G uj x is an ordinary cut locus point, a first conjugate 
point, or both. 

Next we discuss critical distances associated with the boundary ex- 
ponential map, exp 9M , 

exp aM : dM x R — ► M, exp dM (z, s) = 7 Z| „(s), z G dM, 

where v = u(z) is the unit interior normal vector to dM at z. The pair 
(z, s) defines the boundary normal coordinates in M near dM. 
The boundary cut locus distance, r^z), z G dM is given by 

(3) n (z) = max{s > : dist(7^(s), dM) = s}. 
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The set of the corresponding points y = r y z>v {T b (z)) is called the bound- 
ary cut locus, 

UdM = {y e M : y = 7 Z} u(n(z)), z G dM}. 

The boundary cut locus consists of two types of points. We say that 
a point y G ujqm is an ordinary boundary cut locus point if there are 
z, w G dM, z 7^ w with 

r b (z) = T b (w), jz^iniz)) = iw,u{w){n{w)) = y. 

Also, we say that a point y = r y z , v ( z ){T b {z)) G lo x is a focal point if 

the differential, dexp dM \( z>Tb ( z )) '■ T z dM x R — > T y M is degenerate. 
Equivalently, t is a focal point if there is a non-trivial Jacobi field Y(t) 
along 7*,„([0,s]) with F(s) = and F'(0) = W(0), where W is the 
Weingarten map of dM at z. For z G <9M, we define the focal distance, 
Tf(z) to be 

Tf(z) = inf{s > : (iexp^ |( Z)S ) is not one-to-one}. 

Note that y G loqm is an ordinary boundary cut locus point, a first fo- 
cal point, or both. Also, the functions tr, t c , t^, and 7/ are continuous, 
e.g. [261- 

Comparing Jacobi fields Y(s) along the geodesic J z ,u([0, s]) with the 
Dirichlet condition Y(0) = and the Robin condition Y'(0) = WY(0), 
we see that Tf(z) < r c {z,v). Due to the compactness of dM there is 
c > such that 

T c (z, v) > T f( z ) + Co, z G <9M. 

In a similar manner, we can show that tr{z,v) > n(z), z G dM. 
Indeed, assume the opposite, i.e., t = r R (z, v) < r b (z) for some z G dM. 
Denote (y,r)) = (i Z)U {t), -j' Z;U (t)). By duality, T R (y,T}) = r R {z,v) = t. 
Let e > and x £ = j Z;U (—e) = 7 s/fl (t + s). Then 

distjg (x e , y) < t + e < r b (z) + e 

and there is r\ e G S Xe M with y = 7a; e)J)e (dist^(x e , y)). Denote by t E > 
the last time when 7 X£i ^ e (s) hits dM. If e is sufficiently small, we see 
by the short-cut arguments that dist(y,<9M) < T b (z). This contradicts 
the definition of r;, in (3|. 

Due to the compactness of dM, by making Co > smaller if neces- 
sary, 

(4) r R (z, v) > n{z) + c , z G dM. 

Later we will consider intersections of various geodesies on M . In 
these considerations we would like to avoid pathological cases that may 
happen to long geodesies. The first case we analyze is a self-intersection 
of a geodesic. 
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Figure 2. Left: Self-intersection of a normal geodesic. 
Right: Geodesies corresponding to focusing directions. 

Lemma 2.1. Let j z>v , z G dM be the normal geodesic and 

7*,i/(s+) =7*,f(*-)> s + > s_, 
that is, 7 Z)jy intersects itself. Then s + + s_ > 2tr(2:, z/). 
Proof. Assume that 

(5) + <2r fl (z,i/). 

Then s_ < tr(z, z/). Let A = j z>v (s-.), B = 7 2> z,(tr(,2, v)) be points on 
7 Z)1/J see Fig. 2, and denote by Iba = s+ — Tr(,z, z/) the length of the 
"long" geodesic 7«,j/([tr(z, v), s + }). Then, using definition (J2J) of tr, 
s_ = dist(z,A), T R (z,v) — s_ = dist(A, £>), so that the length of the 
broken geodesic 7 Z ,„([0, s+]) U 7«,i/([0, s_]) from z to 2; is 

s + + s_ = dist(s, A) + dist(A, 5) + / BA + dist(A z). 

Since t KjI/ ([s_, Tr(z, v)\) is the unique minimal geodesic between its 
endpoints, Iba > dist(A,B) = tr(z,u) — s_. Therefore, 

s+ + s_ > s_ + {t r (z, v) - s_) + {t r (z, v) - s_) + s_ = 2tr(2, z/), 

which contradicts ((5|). □ 

In the sequel, dists is the Sasakian distance on, depending on the 
context, TM or SM, see [36] . 

Lemma 2.2. Zet e > 0, z e 9M. TTiere is 5 = 5(e) > sizc/z i/ 

(^1,6)^24(^2,6), ie. 7*i,ft(*i) =7« I ,6(*2), h + t 2 = 2t, 
with t < tr(z, v) + 8 and distg((zi, 6), (z, v)) < 8, % = 1, 2 i/jen 



Note that the constant 5 does not depend on 2 6 9M. 
Proof. Assume the opposite, i.e., an existence of points z k e dM, 
(zf, £f) G = 1, 2, z = 1, 2, . . . and a parameter e > 0, such that 



t — ti\ < £ 



1 = 



1,2. 



lim dtot ff ((** (z k ,u k )) = 0, 
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14,^1) = 7** g(*2)> 4 + t h 2 = 2t\ limsup(t fc - r R {z\ v k )) < 0, 

fe^oo 

with t\ — 1\ > 2e. Using continuity arguments and compactness of dM 

we have that there is a subsequence k(p) with 

t; (p) -,r,and 

7*,»/(f+) = 7*,„(t_), t+ + f_ < 2tr(z, i/), t+ - t- > 2e, 

which contradicts Lemma r2.ll □ 

Next we introduce auxiliary functions ^i(z), ^(z), and tm(z), z G 
dM with Hi{z) and // 2 (z) to be determined from the broken scattering 
relation. The function ^i(z) tells when a normal geodesies sent from 
z G M exits M. By the definition of the broken scattering relation, 
R, a point (z, £) G 0+ is in relation with itself, (z, ^)R t {z, £), if and 
only if the geodesic 7^5 ((0, t/2]) on M lies in M int . This makes it 
possible to determine, for any 7^, (z,£) G Q + , its arclength to the 
first hitting point to dM. We denote this arclength by f/,i(z,£) and 
/ii(z) = jii(z, v). 

The function /i 2 (z) is an approximation to Tf(z). If we want to 
determine Tf(z) we can argue as follows: assume that s > Tf(z). Then 
the normal geodesic 7 Z)I /([0, s]) is no longer a shortest path from 7s,„(s) 
to <9M and there are sequences z n — > z, z n 7^ z, s n — > 17 (z), £ n — > r/(z) 
such that 

7^ («n) = 7;w*(*»)> =K z n)- 

In terms of the relation i?, these imply that 

(6) (z, v) R Tn (z n , u n ), T n = t n + s n , 

with s n — > Tf(z), t n — > Tf(z), z n — > z, when n — > 00. 

Therefore, it makes sense to try to find r/(z) using (jHJ). However, 
there are two obstacles. First, it may happen that r/(z) > f-ii(z). 
Second, having (JO]) with z n — > z, T n — > 2i, we want to conclude that 
s n — >■ £, t n — > t. To do so, we intend to use Lemma l2~2l which requires 
t < tr(z, v) which is not known. To avoid these difficulties, we will not 
determine Tf(z) but another function /x 2 (z) that is closely related to it. 

Definition 2.3. Consider the set S(z) of those s G (0, Hi{z)) for which 
there are sequences z n — > z, z n G dM z n 7^ z, T n — > 2s such that 

(7) (z n ,v n )R Tn (z,p). 

Define // 2 (z) = 'vafS(z), if S(z) 7^ and // 2 (z) = /-*i(z) otherwise. 

Observe that /i 2 may be found from the broken scattering relation. 
Lemma 2.4. Function /i 2 : dM — > R + satisfies 

(8) min(//i(z), 7/(z), t r (z, v)) < fx 2 (z) < min(/ii(z), r/(z)). 
and 7,(z) < /i 2 (z). 
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Proof. The right inequality in (jSJ) follows from Definition 12.31 and 
considerations before it. 

To prove the left inequality of (jHJ), let us assume that there is s < 
min(Tf(z), H\(z), tr(z,u)) which satisfies jTj). By Lemma l2~2l applica- 
ble due to T n < 2tr(z, v) for large n, we have 

As s < Tf(z), exp 9M is a local diffeomorphism near (z,s), which con- 
tradicts j9]). This proves (JHj). 

Using definitions /xi and r/, we see by using (J4j) that 

T b (z) < min(-/i 1 (z), r/(z), r fl («, 

This yields r&(z) < ^{z). □ 

Finally, we need a function tm(-z) with tm(z) > Tb(z) having the 
property that, for t < tm(z) the geodesies sent back from a point 
x = J z ,v{t) hit the boundary dM near z in a regular way. Namely, we 
define 

t m (z) = min (jii(z), r fl (2, ^(^))), 2 G <9M. 
As ^(z) < we see by (jlj) that 7^(z) < r«(^). 

2.2. Family of intersecting geodesies. In this section we intend to 
use the broken scattering relation to verify if a given family of geodesies 
intersect at one point. 

Let z G dM, u = v{z ), and x = 'j zo>uo (t ), < t < t m (z ). 
Denote r] = — j' ZOVQ (to). Clearly, r] is the direction of the reverse 
geodesic, ■J Xo ,r]o from xq to z . By considering Jacobi fields along this 
geodesic, we see that the exponential map, exp^ : S Xo M x R + — > M, 
is a local diffeomorphism near (770, to)- 

As t < t r (x ,i]o) and 7so,*?o(*o) hits dM normally, all geodesies j X0)V 
hit dM transversally for r\ G S Xo M close to r] . They determine smooth 
functions z(r)), t{rj) such that ^xa^itijf)) = zijf) G dM. Inverting these 
functions and using transversality, we obtain, in a neighborhood U C 
dM of z a smooth section £(z) : U — > SU and a function t(z) such 
that 

(10) 7z,«*) (*(*)) =X , ZEU. 

In the following, our aim is to determine, using the broken scattering 
relation R, whether, for a given triple {£/,£(•),£(•)} of a neighborhood 
U C dM and functions £(z) and t(z), there exists a point xo G M such 
that 7*,£(z)(£(z)) = Xo for all z E U . 

To this end, we notice that property (fTOl ) implies 

(11) (z,^))^rw (20,^0), M*))^.^^^')), G?7, 
r(z) = t{z) + t , *') = t{z) + t(z'), 
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for smooth £(z), t(z). In addition, 

(12) t{z )=t , dt(z)\ zo = 0, Z(z )=u(z ), 

where the last properties follow from the fact that 7x ,r;o is normal to 
dM. Here, dt(z) = d z t(z) is the differential of the function t : U — > R. 
These observations motivate the following definition: 

Definition 2.5. Let z Q G dM and to > 0. Consider a family T '(zo, to) = 
{U,£(-),t(-)} where U C dM is a neighborhood of z , £ : U -> SM 
«s a smooth section, and t : U — > R is a smooth function. We say that 
jF(z ,t ) ?s a family of focusing directions if £(z), t(z) satisfy condi- 
tions fTJ\) and ( TfljJ , VFe t/ien sa?/ t/iat t/ie geodesies r y z ,e,{z)i z £ U are 
the geodesies corresponding to family J-(zo,to). 

Note that the broken scattering relation R determines if given U, 
£(z), an d t(z) form a family of focusing directions. Our principal tech- 
nical result in this section shows that the geodesies corresponding to a 
family of focusing directions intersect at a single point. 

Theorem 2.6. Let z G dM, t Q < t m (z q ), and J r (z ,t ) be a family 
of focusing directions. Then there is a neighborhood U C U of z such 
that 

lz,t(z)(t(z)) = lz ih v (t ), for all zeU. 

Proof. The proof of this result is rather long and will consist of several 
steps and auxiliary lemmata. 

Step 1. We start with an observation that (fTTTl implies that, for any 
z G U, there are s(z),Hs(z) > such that 

x(z) = Jz,t(z)(s(z)) = Jz ,u (s(z)), s(z) + s(z) =T(z). 

As to < t r(zo, vo), by Lemma [2721 s(z) — > to, 's(z) — > to when z — > z 
and 

(13) s(z ) = s(z ) = to- 
Next we show that s(z),'s(z) are C°°-smooth near z and 

(14) ds(z)\ Z0 =ds{z)\ Z0 = 0. 
To this end, consider the function H(s,z), 

H(s, z) = dist(7 20)I/0 (s), z) + s- T(z), (s, z) G (t - <5, t + 5) x U. 

As to < T^(z , vq), the function i?(s, z) is C°°-smooth a neighborhood 
of (t , Zq) and 

if (t , Zq) = 0, d s H(t , Zq) = 9 s dist(7 20 ^ (s), z )\ to + 1=2. 

Making U smaller if necessary, the equation H(s, z) = has a unique 
solution s = 7s(z) which is C°°— smooth in U with 's(zo) = to- As also 
s = s(z) solves H(s,z) = 0, we see that s{z) = 's(z), z G U. It then 
follows that s{z) = T{z) - s(z) G C°°{U). 
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Let us differentiate the identity H(s(z),z) = with respect to z at 
z = Zq. Due to (fl2l) and the fact that 7z 0ii , is normal to <9M, 

= d z H{s{z),z)\ Z0 = d z s\ zo ■ (9 s dist(7 Z0)J , (s), z )\ s=to + 1) = 2d z s\ zo . 

Thus, <i 2 ?| 20 = and also d z s\ zo = d z (T(z) — s(z))| Zo = 0. 
S^ep 2. Consider the map E G C°°(U; SM), 

E(z) = (x(z),r)(z)) := (j z4(z) (s(z)), -J z ^ {z) (s(z))) , E(z ) = (x ,r)o)- 
Lemma 2.7. The map dE\ Z0 : T ZQ dM — > T X0)m SM has the form 

(15) dE\ zo {v) = (0,ev), veT Zo dM, 

where we identify T xom SM « T X0 M x T V0 (S X0 M) . Furthermore, O : 
T Zo dM — » T^S^M) is bijective. 

Proof of Lemma 12. 7L As x(^) = 7 ZOi ^ (s(^)), it follows from (THl) 
that ofe| 2o = 0, i.e., dE\ Zo is of form ( fT5l) . To show that is bijective, 
observe that 

(16) exp x{z) (s(z)rj(z)) = z, z G U. 

Let us denote Exp(x, £) = exp^, £, (x, £) G TM. By differentiating both 
sides of (fl6l) with respect to z and using cfe| 2o = 0, we obtain 



Step 3. Our further considerations are based on the analysis of the 
intersection of a single geodesic and the geodesies corresponding to a 
family of focusing directions. 

Lemma 2.8. Let z G dM and T{zq, t ) = {U, £( • ), t( ■ )}, t < r M (z ) 
be a family of focusing directions. Let 7(7") be another geodesic in M 
which intersects 7 Z0)i / , 



(17) 7 (0) (n>), 7'(0) + ±7 20 ,^ (r ), r < t m (z ). 



Assume, in addition, that all geodesies "fz,£(z) corresponding to JF(z ,t ) 
intersect 7 near yo, i.e., 



where < r(z) < r\ < tm{zq) and \r(z)\ <i\< inj(M). Then r = to- 
Proof of Lemma 12. 8L Denote yo = 7z ,i/ (?"o)- First we show that 
r(z) is continuous at z . If this is not true, there would be another 
intersection of 7 20jJ , and 7, 




□ 



(18) 



lz£(z)(r(z)) = j(t(z)) 



lz ,u {r') = l{r'), r' < n, r' ^ r , \r'\ < inj (M). 
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This leads to a contradiction as both 7([0, r']) and Jzo,v ([ r o, r ']) are 
unique minimal geodesies between their endpoints. Thus r(z) is con- 
tinuous at Zq. 

To prove the claim, we assume that r 7^ t . Our next goal is to show 
that the map * : [/ x 1 + -> M, 

#(*,r) = exp z (r£(z)) 

is a local diffeomorphism near (z ,ro), see the right part of Fig. 3. 
Indeed, as t ,r < Tr(zq,Vo), the map exp is a local diffeomorphism 
near (t - r )r) , where x = 7»,n>(*o)> Vo = (*(>)■ Thus, 

rfex Px \(to-r )vo '■ T (to-ro)rio( T xo M ) T S/o M 

is bijective. Using the definitions for s(z), = (x(z),T)(z)) intro- 

duced earlier we have 

(z, r) = 7 B (a)(s(z) - r) = exp x(z) ((s(z) - r)^(^)). 

By (EED and (fT3|) . cis(z)| Zo = and s(z ) = to, which together with (IBj) 
imply that 

^l(^o,r )(C,P) = dexp^ Icto-roJiwCCto - ^o)©C - PPo) 

for ( G T ZQ dM and p G M. Thus, by Lemma 12.71 and bijectivity of 

^ ex Px I (to— nj)%! 

d*| (zo , ro) : T Zo dM xl-» T^M 

is bijective, i.e., \I> is a local diffeomorphism near (z , ^o)- 

Now, let S be an (n — 1)— dimensional submanifold which contains 
a part 7(— e, e) of 7 near y an d is transversal to •y zo ,u at y , see Fig. 3, 
the existence of such submanifold guaranteed by (TT71) . Introducing the 
boundary normal coordinates (w,n) associated to S, with n = on S, 
we rewrite \I/ in these coordinates as 

^(^r) = (w(z,r), n(z,r)). 

By transversality, ^(z ,r ) 7^ 0. This implies that for any z near z 
the equation n(z, r) = for has a unique solution r = r(z). Moreover, 
r(zo) = r and the function r(z) is smooth in a neighborhood of z$. 

Now r(z) and r(z) are continuous at zq and they both solve the 
equation n(z,r) = 0. Thus, there is a neighborhood U C U of z such 
that r(z) = r(z) for z E U. As also \I/ is a local diffeomorphism, we 
see that if U is small enough, then ^ : U — > ^(C/) C S, where ^(2) = 
\P(z, r(;z)), is a diffeomorphism of (n — l)-dimensional submanifolds. 
On the other hand, condition (0781 ) implies that *&(U) C 7(— s,e). As 
7(— e, e) is a one-dimensional submanifold of £, we get a contradiction 
for n > 3. Thus, r = to- n 

Step ^. Let < e < imin(inj (M),r R (z 0l u) - t ) and < 5 < 5(e) 
where 5(e) is defined in Lemma [2721 We choose a neighborhood U C U 
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Figure 3. Left: Submanifold £ contains geodesic 7 and 
is transversal to j ZOt y. Right: Geodesies corresponding to 
^(zQjto) almost intersect at the point xq = t^^o) and 
define coordinates near p = 7 20i ^(r ). 



of z so that 

\t(z)—to\<e and ds((z,£(z)),(zo,Uo)) < 6 for z G U. 

By Definition 12.51 there exist functions si(z, z'), s 2 (z', z) > 0, z, z' G 
U, such that 

7*,C(«)(si(z,^)) = lz>,t(z>)(s2(z',z)), s x {z,z') + s 2 (z',z) = t(z) +t(z'). 

By Lemma [2721 these imply that 

(19) \t - Sl (z,z')\ < 2e, |to-s 2 (/,z)| <2e. 

Consider a geodesic 7(5) = j z ',£(z'){ s + ^(^'j ^o)) for some fixed z' G U , 
z' 7^ 2; - It follows from ( fT9"l) that Lemma 12.81 is applicable to the 
family jF(2 ,to) an d the geodesic 7 with r x = r R (z , z/ ) — 2e, i x = 2e. 
Thus, 7^', S ( Z ') and 7^ intersect at x = 7^0(^0)- As z' e U \ {z } is 
arbitrary, all geodesies corresponding to family T(zq, t ) with a starting 
point z' G U intersect in x . □ 

Later on we will need the following modification of Lemma [2751 which 
do not require that all geodesies of JF^q, t ) intersect 7 near y . 

Lemma 2.9. Let z G dM and T{z§, t ) = {U, £( • ), t( ■ )}, t < r M (z ) 
be a family of focusing directions. Let 7(7-) be another geodesic in M 
which intersects all geodesies "fz,£(z) corresponding to JF(,2 , t ) ; 

7*,«*)(r(*)) = 7(t(z)), 

where < r(z) < ri < r^zo) and |r(z)| < L, where L > is arbitrary. 
Assume, in addition, that h(z) = r(z) + t(z) is continuous. Then 
Jz,£(z)(t(z)) = j(h(zo) — to) when z is sufficiently close to zq, i.e., all 
geodesies intersect at the same point. 



14 



YAROSLAV KURYLEV, MATTI LASSAS, AND GUNTHER UHLMANN 



Proof. We first show that there is only a finite number of intersec- 
tions of 7z ,i/ ((0, ri)) with j([—L,L\). Let 7i, . . . , Tjy G [—L,L] and 
Tq, . . . ,Tq G (0, ri) define the points of the intersection, 

As all geodesies in balls of radius inj (M) are shortest and r ] < r 1 with 
lz ,u ([0, r i}) being the shortest between its endpoints, 



N < 



2L 



inj(M)J 



+ 1, 



where [t] denotes the integer part of t G R. 

Let < e < |inj (M) and U(p) = dMf]B(z , p), where B(z , p) C M 
is the ball with center z an d radius p. Then there is p > such that 

min \r(z) — rl\ < e, for z G U(p ). 

Indeed, otherwise there is a sequence z n — > zo with r(z n ) — >■ r < 
tr(2;o, ^(^o)) and r(z n ) — > r, |r| < L, such that 

7^,^o (r) = 7(t% i = !> • • • > N, 

which is a contradiction. 
For < p < po, denote 

Vj{p) = {z eU{p) : 7*,C(*)W = 7(t), r (^) + T ( z ) = h ( z ), \ r ~ r o\ < £ }- 

Sets Vj(p) are relatively closed U(p) and, therefore, measurable on 
dM. As \Jj =1 Vj(p) = U(p), we see that for some j the set Vj(p) has 
non-zero (n — l)-dimensional measure. However, if ^ to, the same 
considerations as in the proof of Lemma 12.81 by replacing r by r 3 Q and 
using a relatively open neighborhood U C Vj(p) of show that the set 
Vj(p) has (n — 1)— dimensional measure equal to when p > is small 
enough. This shows that there are j and p > such that r J = t and 
\ Vj(p) has (n — 1)— dimensional measure equal to 0. Thus Vj(p) 
is dense in U(p). As e > is arbitrary, the continuity of the geodesic 
flow shows that 7^ ,^ (to) = 7(^(^0) — ^o)- Together with Theorem 12.61 
this completes the proof. □ 

In the following we say that two geodesies pit) and pit) coincide if 
p{ti) = p(t 2 ) and p'{ti) = ±p'(t 2 ) for some ti,t 2 G R. Note that this is 
equivalent to p{t) = p(a + t) or p{t) = p(a — t) for all t in a non-empty 
open interval and a G R. 

2.3. Reconstruction of the boundary cut locus distance. 

Lemma 2.10. The boundary, dM, and the broken scattering relation, 
R, determine the boundary cut locus distance Tj,(z), z G dM. 



RIGIDITY OF BROKEN GEODESIC FLOW 



15 



Proof. We recall that for t < r b (z ) the point z in the unique point 
of dM closest to x = ■y zo ,v (t ). On the contrary, when t > T b (z ) 
there is another point w G dM with dist(7 ZOjl , (to), w) < to- What is 
more, considerations in the beginning of Section 2.2 show the existence 
of a family T[zo, to) of focusing directions for t < tm{zo). Recall that 

n(z ) < t m (z ). 

Thus, when T b (zo) < t < Tm(zo), there is a family jF(2 ,t ) = 
{U, £(■ ), t(- )} of focusing directions, a point w G dM, w ^ Zq, and 
s < t such that 

(20) {z, R t(z)+S0 (w, u(w)), zeU. 

Our next aim is to show that when t < T b(zo), there are no w G dM 

and jF(zo,to) satisfying (1201) with so < to- 
Assuming the opposite, there is a neighborhood U C dM of zq and 

a function r(z) with 

(21) 7/z£{z){r{z)) = J w , v (w){t(z) -r(z) + s ), z G U. 
Next we prove that 

(22) r = limsupr(^) < t . 

Assume that (1221) is not true. Then there is a sequence z n — > z with 
r(z n ) — ■> r > to- By the continuity of the exponential map, it follows 
from ([21]) that 7« ,^( r o) = lw,u(w)(t - r + s ). Thus, by the triangle 
inequality, 

dist (w, 7z ,v {to)) 

< dist(w,^ WjV(w) (t - r + s )) + dist(7 ZOji , (r ), 7* 0) fo(*o)) 

< (*o - r + s ) + (r - t ) < s < t , 

which contradicts the definition j3|) of r b . Thus (l22l) is valid. 
Therefore, by making U smaller if necessary, we have 

r(z) < T M (z ), z eU. 

Assume first that geodesies 7 20i y and jw^w) do not coincide. Applying 
Lemma l2~9l with j(t) = y w ,v(w)(to + s o — r o + T ) and L = 2t , we obtain 
7^0,^0(^0) = Jw,v{w)( s o)- As s < to this contradicts with the definition 

Of T b . If 

"f zo ,vo an( l Iw, v(w) coincide, condition w 7^ z$ implies that 
7^0,^(20) (^0 + so) = w. Then we would have dist(xo, dM) < dist(xo, w) < 
s < Tb(zo), that is not possible. 

Finally, by Lemma EH the relation R determines the function 
satisfying r b (z) < ^{z). Let J(zo) be the set of those to G [0, /^(-^o)] 
for which there are w G dM, s < t , and jF(z ,t ) satisfying f[20l . 
If r b (z ) < fi2( z o), we see that (r b (z ), ^(zo)) C J(z ). Thus we can 
determine r b (zo) by setting T b (z ) = inf J(z ) if J(z ) 7^ and T b (z ) = 
jj.2i.Z0) otherwise. □ 
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2.4. Boundary distance representation of (M,g). Next we con- 
struct of isometry type of manifold (M, g) by showing that the broken 
scattering relation, R, determines the boundary distance representation 
H{M) of (M, g) that is the set 

K{M) = {r x : xe M} c C{dM), 

where r x : dM — > R are the boundary distance functions 

r x (z) = dist(x, z), z G dM. 

It is well-known, e.g. EH EE] that the set 1Z(M) possesses a natural 
structure of a Riemannian manifold with the map 

TZ:M^TZ(M), 7Z(x)=r x (-) } 

being an isomorphism. What is more, this metric structure can be 
identified just from the knowledge of the set 1Z(M). An additional ad- 
vantage of dealing with TZ(M) is the existence of a stable procedure to 
construct a metric approximation, in the Gromov-Hausdorff topology, 
to (M, g) given an approximation to 1Z(M) in the Hausdorff topology 
on L°°(<9M), [23j. To construct 1Z(M), we assume that the function 
U is already known. We start with finding dist^M on dM which is 
inherited from (M,g). We define that distgM^i, z 2 ) = oo when z\ and 
z 2 lie on different components of dM. 

Lemma 2.11. The boundary, dM, and the broken scattering relation, 
R, determine, for any Zi,z 2 G dM, the distance distgM( z i, z 2) along 
dM. 

Proof. It is enough to consider the case when z\ and z 2 are in the 
same component of dM. 

Using boundary normal coordinates, we see that there is e > and 
c > such that 

(23) | dist (2/1,2/2) - dist 9M {yi,y2)\ < c e 3/2 , 

if distal (2/1,2/2) < £ 3/4 , £ < £o- Let x 2 = J y2 ,v 2 (e 5/4 ) . Making e > 
smaller if necessary, we see that there is a unique shortest geodesic in 
M, 7 ai ,f 15 with (2/1, £1) G fLf, from yi to x 2 . Moreover, using again 
boundary normal coordinates, we see that 

(24) |dist(t/i,x 2 ) + dist(x 2 ,j/ 2 ) ~ dist dM (y 1 ,y 2 )\ < Ci£ 5/4 . 

Let /i = /i([0,/]) be a shortest geodesic of dM from z\ to z 2 . Let 
N eZ+,e = l/N and y d = n(ej),j —0,...,N. Define Xj = 7 %) ^(e 5/4 ) 
and associate with each j = 1, . . . , N a broken geodesic atj which is the 
union of the geodesic from to Xj and from Xj to j/j. Inequality ([241) 
implies that if — >• 00, then 

AT 

(25) |dist 9M (^i, ^2) - ^ (dist(j/j_i, Xj) + dist(y i; Xj)) \ < c 2 e 1/4 -> 0, 
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Motivated by this, define for iV G Z + and e = 1/N 

N 

d N {z 1 ,z 2 ) = inf ^Sj, 

3=1 

where the infimum is taken over the points yj G dM, j = 0, 1, . . . , N, yo = 
Zi,JJn = z 2, which satisfy the following condition: For any j = 0, . . . , N— 
1, there are r/j G S yj M, (vj,r)j) g > and positive Sj < e 3//4 such that 

((Vj, %)> (Vj+i, K%'+i))> e R, j = 0, 1, . . . , N - 1. 

Using fl23l) we see that £^(21,22) > dist^C^i, z 2 ) — Cse: 1 / 2 . On the 
other hand, as we saw in f l25l ). there are yj,T]j, and Sj such that 

Idist^Af {zi, Z2) — cIn{zi, z 2 )\ < C4£ 1//4 = ciV~ 1//4 — >• 0, when iV — > 00. 
Thus we get that 

dM, dM {zi,z 2 ) = lim d N (z 1 ,z 2 ). 

N—*oo 

□ 

Next we determine the distance between boundary points with re- 
spect to the metric g in M. 

Lemma 2.12. The boundary, dM, and the broken scattering relation, 
R, determine the distance function dist(xi,x 2 ) for x\,x 2 G dM 

Proof. By [T], for any X\, x 2 G dM a shortest path connecting them is 
a C 1 — path. Let x(s), s G [0,1], I — dist(a; 1 , x 2 ), x(0) = X\, x(l) = x 2 be 
such a shortest path, parameterized by the arclength, that connects x\ 
to x 2 in M. Moreover, by [TJ it holds that if x(s) G M mt for s G (a, b), 
then x((a, b)) is a shortest geodesic between x(a) and x(b) in M. 

Clearly, the set of s G [0, 1} such that x(s) G M int is open. By (j25j>. 
for any e > there is a finite number points Oj, i = 1, . . . ,p, a p+ i = /, 
and bi, i = 1, . . . ,p with < ai < 61 < a 2 ■ ■ ■ < b p < a p+ \ = I such that 
Zi = x(ai),yi = x{pi) G dM and 

(26) dist(xi, x 2 ) < distal (^i> z\) + 

+ {^L dist (^ 5 Vi) + dist dM (yi, z i+ i) \ < dist(xi, x 2 ) + e 

and there are shortest paths 7« i) » 7i ([0, in M of length = 6j — a* 
from to j/j that satisfy 7z ij77j ((0, 6j — a,)) C M mt . Next we will relate 
(l26l) to the broken geodesic relation. Recall that relation R involved 
broken geodesies that start and end non-tangentially to the boundary. 
Because of this, we consider for tangential 77, the vector & = (1 — 
h)^ 2 r]i + h l l 2 i>(zi) G S Zj M. If 77^ is non-tangential, we set & = 77^. 
When /i > is small enough and < U is sufficiently close to k, 
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we have that 7 Zi) ^((0, Sj]) C M mt , and the closest boundary point to 
7zi,&( s i)j denoted iji, satisfies 

dist (7^(si),2/i) < ^ dist 9M (yi, yi) < J. 

Consider the broken geodesic from 2$ to ^ which is the union of the 
geodesic from Zj to 7 Zi ,f i (si) and from 7^,^(sj) to ifj. It has the length 
ti < k + s/p and non-tangential starting and ending directions. Thus 
{zi^^R^iyi^f). These considerations show that 

I P 

dist(xi, x 2 ) = inf I dist^M^i, z 1 ) + ( s ^t i + dist dM {yi, z i+1 )) 

\ i=i 
where the infimum is taken over ti > 0, Zi,y~i G dM, and directions 
£i,C« such that z p+ \ = x 2 and the relations (z iy £i)R tl (yi, Q) are valid. 
□ 

Theorem 2.13. The boundary, dM, and the broken scattering rela- 
tion, R, determine the set 1Z(M) C C(dM). 

Proof. Let ujqm be the boundary cut locus on M. As M\ujg M is dense 
in M, it is sufficient to find TZ(M\uqm)- Recall that, for x G M\ujqm, 
we have x = 7 ZOiJ / (£o)i where t = dist(x ,<9M) < r b (z ) and z is 
the unique boundary point closest to x . Using the broken scattering 
relation R, we intend to determine, for any Wo G dM, D(z ,t ,Wo) '■— 
dist(xo, wo). 

Let x(s) be a shortest path from xo to wo parametrized by the ar- 
clength. Denote by w = x(sq) the first point where x(s) is in dM. 
Clearly, 

(27) dist(a;o, w ) = s + dist(u>, w ), s > t . 

By [TJ, the path x([0,s ]) is a geodesic in M. We denote r\ = —x'(s ) 
so that x = 7 Wi?? (so). As t < r b (z ) < t m (z ), there is a family of 
focusing directions ^(zo^o) = {U, £( -),t(- )} such that for Si = s , 
Wi = w, and 7/1 = r\ we have 

(28) (w 1 ,r ]l )R Sl+tiz) (z,^(z)), z G U. 
After these preparations we will show that 

(29) D(z , t , w ) = inf (dist(tu , Wi) + «i) 

where infimum is taken over W\ G dM, 771 G S Wl M, and Si > t such 
that there is a focusing sequence JF(z ,t ) = {£/,£( • ),t( ■ )} satisfying 

(EH}. 

Formula ( |27l) shows that the infimum on the right side of (1291 is less 
or equal to D(zo, to, wo)- Thus to prove (i29l . it is enough to show that 
if wi,rji, and si satisfy (i28l then p = dist(u>o, u^) + s\ > dist(xo, wo). 

Assume now that ([28]) is valid. Then, for some r(z),r(z), r(z) + 
t(z) = si + t{z), we have that j z ^ (z) (r(z)) = 7(r(z)). 
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Keeping aside the trivial case when the geodesies 7 ZOji , and 7 M1)J?1 
coincide, consider first the case when limsupr(z) — r > t . Denoting 
7z ,^o( r ) = x i5 we then have 

dist^i, x ) < dist(wi, Xi) + dist(^i, x ) 

< (si +t - r) + (r - t Q ) < st, 

yielding p > dist(w ,Wi) + dist(wi, x Q ) > dist(w ,x ). If, however, 
limsup 2 ^ Zo r(z) = r < t 0l we are in the situation of Lemma [2T9l which 
shows that 

7*d,h>(*o) =T(«i)> 

yielding again that p > dist(w , x ). □ 

As the set TZ(M) can be naturally endowed with a differential struc- 
ture and a Riemannian metric so that is becomes isometric to (M,g), 
see e.g. [24, 25], we have finished the proof of Theorem 11.11 □ 



3. Proofs for the radiative transfer equation. 

3.1. Notations. Let X be a manifold with dimension n and A 1 C 
T*X\0 be a Lagrangian submanifold. Let ( ) of 

be local coordinates X with x' = (xi, . . . , x^), x" = (x^+i, • • • , Xd 1 +d 2 )i 
x" 1 = (x^+da+i, • • • , x n ), and (f>(x, 9), 9 G M. N be a non-degenerate phase 
function that parametrizes We say that distribution u G V\X) is a 
Lagrangian distribution associated with Ai and denote u G I m (X; Ai), 
if it can can locally be represented as 

u{x)= [ e i ^ e) a(x,9)d9, 

where a(x, 9) G S m+n ' A - N l 2 {X xR N \ 0), see [El EES, EHJ- 

Let Si C I be a submanifold of codimension d\. We denote its 
conormal bundle by N*S = {(x,£) G T*X \ : leS.^l T X .S}. If 
^ = { x ' = 0} in local coordinates, Ai = JV*5i and u G J m (X;Ai), 
then locally 

u{x) = [ e ix ' e a(x, 9') d& , a(x, 9') G S"{X x R dl \ 0) 

where p = m - di/2 + n/4. We denote I m (X; Si) = I"{X; N*Si) and 
say that / M (AT; Si) are the conormal distributions in space X associated 
with submanifold S\. 

Also, we denote by I p,l (X; Ai, A 2 ) the distributions u in £>'(X) asso- 
ciated to two cleanly intersecting Lagrangian manifolds Ai, A 2 C T*X, 
see [Tj2 [29]. Let Si and S 2 be submanifolds of M of codimensions 
di and di + d 2 , S 2 C Si. If in local coordinates Si = {x 1 = 0}, 
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S 2 = {x' = x" = 0}, and Ai = N*S U A 2 = N*S 2 , then the distribution 
u G I p,l (X; Ax, A 2 ) can be locally represented as 



where a(x,9',6") belongs to a product type symbol class S fl ''^"(X x 
(R dl \ 0) x R d2 ) containing symbols a G C°° that satisfy 



d}[d%,4,a(x,6',9")\ < C a p lK (l + \9'\ + + |0"|r H/3 ' 



for all x £ K, multi-indexes a,/3,j, and compact sets K C X. Above, 
/j = p + l - di/2 + n/A and /jf — —I - d 2 /2. 
By [IH [29], microlocally away from Ai fl A 2 , 



P> l {A ,A 1 )cI p+l {A \A 1 ) and J^(A , A x ) c F(Ax \ A ). 



Thus the principal symbol of u G / P ''(A , Ax) is well defined on A \ A 1 
and A x \ A . 

3.2. Born series. In the sequel, we denote the distance on (N,g) by 
d(x,y) = dist(x, y). Let Jx^it) be the geodesic on (N,g) with initial 
point x and initial direction £ G S XQ N. Denote 



The measurement operator A can be extended to distributions w sup- 
ported in SU. In the following we consider u corresponding to w (x, £) = 
8(x ,€q)( x i£)) x o £ U. We assume that 7 XO) ,t (R + ) intersect the strictly 
convex manifold M C N '. To analyze the corresponding solution, let 
us denote the specific geodesic on which the leading order singularities 
propagate by 70 = 1x ^ - Also, we denote the corresponding spray in 
SN by r]o = rj Xo ^ . 

Let Uo(x,£,t) be the solution of the equation (pQ) with S being zero, 
that is, Hu + au = 0, w |*=o = w . Then u (t) = c (x)5 vo{t) (x,£), 
where Cq(x) is a non- vanishing smooth function. To simplify notations, 
we consider the equation for all t G R, obtaining 



In the following we analyze the higher order terms in the Born series, 
that is, 




{j x4 (t) eN: tern}, 

{(T*,e(*), e SN : t G R}, 

{(7x, f (t),7^(t))G5iV: tGR+}. 



u (t,x,£) = co(x)6 Vo (t)(x,£), (t,x,£) G R x SW. 



= QSuj-i, j > 1, 



where Q is defined by v — QF where 



Hv + av = F in R + x SN, 



v\ t =o = 0. 
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We note that there are Ci,C 2 > so the solutions u w of equation 
(Cp) satisfy 

(30) \u w (x,£,t)\ < de^Wwh^sN). 

To analyze the singularities of w, let us take the Laplace transform 
in time t and consider u(k,x,(,) = (u(- , x, £))(&). By (l30l ) the Laplace 
transform is well defined for fc 6 C, Refc > C*2. In the following, we 
consider k first as as a parameter, and denote u(x, £) = «(a;, £, fc). Then 

(k + H)u + cru — Su = w in (x, £) G SW, 
where w (x,£) = <5( xo ,£ )(x, £) and 

^,0 = -e'|^-^(a;)- 9 '" 



The operator H + k : C°°(SN) -> C°°(SN) has Q fc a parametrix, see 
[321122], that satisfies (QF)(k) = Q k (F(k)). Also, we denote = 
Uo(fc) + u sc (k), where u sc {k) = u\{k) +U2{k) + . . . . 

Consider now a Born iteration starting at a general Wo(k). Since the 
coefficients of H are smooth functions and the kernel of S is a smooth 
compactly supported function, we that for any s > there there is 
C-i = C 3 (s) > such that for Rek > C 3 the Born series 

oo 

(31) w(k) = ^(QkSy^woik) 

converges in Sobolev space Hi 0C (SN) when w (k) G H[ oc (SN). 

3.3. Properties of the compositions of the operators S and Q^. 
Lemma 3.1. We can write S = S1S2, 



Sjf(x,C) 



[ Kjix^ftfix^dSg), j = 1,2 



where £, £') G C™(SN x SN). 

Proof. Interpreting 1 as a parameter, we define K x : L 2 (S n ~ 1 ) — > 
L 2 (S n - 1 ) by 

KJ(0= ! K(x,£,Of(OdS(0- 
Js"- 1 

As the kernel K(x, £, £') is smooth, we see that for all a G N n and 
l,m G N there is a constant c a ; m such that 

(32) sup ||^(1- A^ m K(x^,0\\c^>cs^) <c a im, 



n-1 



where is the Laplace-Beltrami operator of the (n — l)-sphere S 
Let a m > be numbers such that < a m < e~ m min(l, c~ lm ) for all 
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max(|a|,/) < m. Then the operator 

oo 

B = 53o ro (l-A e ) rn 

m=0 

defines an unbounded non-negative selfadjoint operator B : L 2 (S' n_1 ) — > 
L 2 (S' n ~ 1 ) having an inverse J = B~ l that can be extended to a smooth- 
ing operator V\S n ~ l ) — > C ac (S n ~ 1 ). Moreover, by ( 1321 ) we see that for 
any x the operator L x = BK X defines a smoothing operator V\S n ~ l ) — ► 
C°°(5' n_1 ) and its Schwartz kernel At (£>£') is a C°°-smooth in all vari- 
ables (£,£,£'). Thus we prove the assertion by defining K 2 (x,£,£') = 
L x (Z,£ r ) and K^x,^,^) = J (£,£'); where J(f,0 is the Schwartz ker- 
nel of J. □ 

The Born series iteration can be written as 

Uj(k) = Q k S 1 A j ~ 1 S 2 u (k) 

where A = S 2 QkSi- To analyze the operator A we consider first the 
case where K(x,£,£') would be the constant 1. Denote by S c the 
operator corresponding to a constant scattering kernel K(x,£,£') = 1. 
For this purpose, we introduce operators T = 7r* : L 2 (SN) — > L 2 (N) 
and T* = 7r* : L 2 (iV) -> L 2 (SiV), that is, 

Tu(x) = c~ 1 u(x,£,)dV g (£), T*v(x,£) = v(x), 
Js x n 

where c n = vol(S' n_1 ) and V g is the volume on S X N. 

Lemma 3.2. Let Z = SN x SN, L = {(x,£,y,r)) G Z : x = y}, and 
S = N*L . The Schwartz kernels of A c and A satisfy 

(33) A c (x,C,y,v) e r\Z-L Q ) = r(Z;S ), 

(34) A(a;,e,l/,r/)G/ p (^So) 

where r = — (n + l)/2, p = r + and £ > 0. 

Proof. Clearly, TT* = I and S c = T*T. Thus we have S c = SfS% 
where S% — S 2 — S. In the local coordinates S c has the Schwartz kernel 

S c (x, £, x', £') = S{x - x') E I\Z- L ) = J mi {Z; E ), 

where m-i = (1 — ra)/2. To analyze A = S 2 Q k Si, we first consider the 
operator 

A c = S c 2 Q k Sl = T*TQ k T*T. 
Denote Q k = TQ k T* : L 2 (N) -> L 2 (iV) and let u e Cg°(iV). Then 

(Q fc T*w)(x,£) = / /i(ar,^,s,A;)T;(7 a . >f (s))d5 
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where h(s, x, £, k) is the solution of the differential equation 

(35) d s h(s, x, £, k) + (k + a(j x> ^(s)))h{s, x, £, k) = 0, 
h(s,x,£,k)\ s=0 = 1. 

Note that 

(36) h(x,Z,s,k) = e- ks h(x,^ Sj 0). 

Thus, using the assumption that the manifold N is simple, we have 

(37) (TQ k T*v)(x) = [ [ h(s,x,Z,k)v( lx ,z(s))dsdV g (t) 

[h(s(x, y),x, £(x, y), k)j(x, y)]v(y) dV g (y), 



N 



where s(x,y) G (— oo,0] and £(x, y) G S X N are defined by exp~ 1 (y) = 
s(x,y)£(x,y), and j(x,y) = det(dexp x {y)" 1 is the Jacobian determi- 
nant where dexp x \ y is the differential of the map exp x evaluated at y. 
Since (N, g) is simple, the kernel b(x, y) := h(s(x, y),x,£,(x, y), k)j(x, y) 
is smooth outside the diagonal and behaves near the diagonal as 

b{x,y) ^ e- kd ^ y) d{x,yf- n . 

Using (1371) we see that Q k is a pseudodifferential operator of order (—1) 
(for a similar argument see [37]). 

The Schwartz kernel Q k (x, x') G I~\N x N; diag (N x TV)) of Q can 
be written as 

Q k (x, x') = [ e i{x - x ' ye a(x, x', 9)d9, a G S~ 1 (N X N x R n \ 0). 



The same expression defines a function Q k (x,£,,x',!;') := Q k (x,x') G 
I-\SNxSN; L Q ). This function is the Schwartz kernel of A c = T*Q k T 
and thus we see that the first part of the assertion, the formula (1331 ) is 
satisfied. 

Next we consider the Schwartz kernel of A, that is, A(x, £, y, rj). It 
can be written as a product 

A{x,i,V,ri) = AC (^,^,y,v)J(x,^y,ri) 
where (using the Riemannian normal coordinates at x) 

\y — x \ \ x — y\ 

Now Ji(x,y,z) := Ki(x,£,z/\z\) and J 2 (x,y,z) := K 2 (x, z/\z\, £) 
are homogeneous functions if degree zero in z, and we see that [T5l 
formula (1.2)] 

K 2 (x, £, f^), K^y, -^L,Tj) G r n (Z; L ). 
\y — x\ \x — y\ 
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Now we can write A as the product of K\, K 2 , and A c . To analyze this 
product, we need the following lemma extending results of [T5] for less 
regular conormal distributions. 

Lemma 3.3. Let Z be a manifold of dimension d and L Q be a sub- 
manifold with codimension n. Assume that A G I~ d (Z; L Q ) and B G 
I^(Z\Lq), fj, < 0. Then the pointwise product AB G /^(Z; L ) for 
any e > 0. 

Proof. Let (V, z") be local coordinates of X such that L = \z! = 0}. 
Then 

A{z)= [ e iz ' e a(z,9)d9, B{z) = [ e iz ' e b(z,9) dO, 

JR d JR d 

where a{z, 9) G S- d (X xR d \0) and b(z,9) G S»(X x R d \ 0). The 
symbol c(z, 6) of the product A(z)B(z) is given by the convolution 

c(z,6) = [ a{z,6-6)b{z,6)dO, 

and a simple computations shows that 

|c(*,0)|<c / {\ + \e -e\Y{\ + \e\)- d de <c\i + \e\Y +£ , 

with e > 0. Indeed, decomposing the domain of integration as M. d = 
B{0, i|0|) U B{6, l\6\) U (R d \ {B{0, ||0|) U B(6, we see that 

\c(z,9)\ < 0^^191+ C 2 \9\- d \9\ d+ ^l + 5^ d \og\9\) + C 3 \9r 
< C'(l+ \9\Y +£ J 

where \9\ > 1 and is one if \i = —d and zero otherwise. The 

derivatives of c(z, 9) can be estimated in similar way, and we obtain 
that c(z, 9) G S^ +£ {X x R d \ 0). ' □ 

Lemma 13731 for the product of Ki, K 2 , and A c implies ( 1341) . This 
proves Lemma [3721 □ 

The previous result says, roughly speaking, that A is like a \l/DO of 
order (—1) operating in (x, y)-variables when ^ and r\ are considered as 
parameters. 

Next we consider powers of A. Next, S' denotes the canonical rela- 
tion corresponding to the Lagrangian manifold S . We see that Sq x Sq 
intersects cleanly T* SN x diag (T* SNxT*SN) x T* SiV with the excess 
d = (n - 1). Thus using [121 Thm VIII.5.2], we see that 

A 2 = A o A G J- 2 ^/ 2 (Z; S ) = I P2 (Z- S ), 

where p 2 = — (n + 3)/2 + 2e with any e > 0. Iterating operator A, we 
see that 

A> eIH(Z;X ) = r 1 - 3+s (Z ] L ), Pj = ±- - j + e, e > 0. 
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3.4. Singularities of the terms in the Born series. In the follow- 
ing, let A = N*Y and A 1 = N*(Yi), where 

= {(io(t),Ut))eSN: teR}, 

Y x = S(-y ) = {(x,Z)eSN: xe l0 {R)}. 

Moreover, let P = P(x, £, D x , D s ) = H + k, 

char (P) = {(x,£,x,0 G T*(SN) : £% + C^U^k = 0}, 

and let B.(x, £, x, £) be the bicharacteristic of P(x,£, D x , (i.e. inte- 
gral curve of the Hamilton vector field in T*(SN) \ 0) starting from 
(x,£,x,£) G T*(SN). Then the flow-out canonical relation generated 
by char (P) is 

A' P = {(x,Z,x,ly,(,y,Oe (T*(SN)\0) x (T*{SN)\0) : 

(x, f , x, f) G char (P), (y, C, y, C) e H(ar, f , J, 0}. 

The flow-out of Ai in char (P) is the Lagrangian manifold A2 C T s N\0 
satisfying A 2 = A' p o A' . 

Lemma 3.4. We have 

Uo(x,Z,k) = Co{x,k)S Vo (x,C) G F°(SN;A ), 

where Co(x,k) is a smooth non-vanishing function andr = (2n — 3)/4. 
Por j > 1, 

1 72 

(38) Uj{k)e n'S{SN) A U A 2 ), r j = ---j + e6 j > 2 , e>0, 
where 5j> 2 is one if j > 2 and zero otherwise. 

Proof. For the zeroth term in the Born series the claim is true by 
definition. Next we analyze the higher order terms. Clearly, 

S 2 u (x,£,k) = K 2 (x,£,r)(x))(S c u )(x,£,k), 

where r)(x) G S X N defines a smooth vector field such that if x = 7q(s) 
then r)(x) = jo(s). A simple computation shows that A' x S' intersects 
diag{T*SN x T*SN) x (T*SN) transversally. Now S 2 G I°{SN x 
SN; L ) = I mi (SN x SN; S ), where m t = (l-n)/2 and by pi Thm 
25.2.3] that S 2 can be considered as a continuous operator 

S 2 :F°(SN;A )^I S (SN;A 1 ), 

where s = r + mi and A[ = AqoSq. A simple computation shows that 
A^oE^ = A' 1? and that A[ x £' intersects diag(T* SNxT*SN) x (T*SN) 
cleanly with excess e = (n — 1). Thus we have by [HI Thm 25.2.3] that 

A j S 2 u {k) G J« +mi+e / 2 ( i S'A^;A 1 ). 

Again, as A[ oT,' = A[, and Ai x E' intersects diag(T*S7V x T*SN) x 
(T*SN) cleanly with excess e, we see that since Si G I mi (Z; E ), 

(39) S 1 A j S 2 u {k) G pj+^i+e/2)^ SN . Ai ) = jp 3 ( SN . Al ). 
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To analyze Uj(k) = Q k S\A^~ l S 2 u (k) , we observe that the operator 
P = H + ik is a first order operator of real principal type. As Qk is its 
parametrix, it follows from [29] that the Schwartz kernel 

(40) g fc e/i- 1 '-?(Z;A r . z ,A P ), 

where A' T « Z is the diagonal of T*Z x T*Z and A' P C T*(Z) is the flow- 
out canonical relation generated by char(P). Now N*Yi intersects 
char (P) transversally. Hence we obtain (1381) by [15], Prop. 2.1]. □ 

3.5. Principal symbol of the singularity. For any s > there is j 
such that Uj (k) G H[ oc (SN). Using the convergence of the Born series 
(13T1) . we see that the series Uj (k) +Uj 0+ i(k) + Uj 0+2 (k) + . . . converges 

Next we consider how to find the geodesic 70 in U. To this end 
we observe using (1391 ) that Tu{k) = Tuo(k) + Tu sc {k) G J°(iV;7o) 
and Tu (k) G I°(N; 70) have the same non- vanishing principal symbol. 
Thus Tu{k) in Z7 determines U D 70- 

Moreover, the above convergence of the Born series in Sobolev spaces 
and ([381) yield that «i(A;) and u sc {k) = Ux(k) + u 2 {k) + . . . are both 
elements in I ri ~i(SN] Ai, A 2 ) and they have the same principal symbol 
on A 2 \ Ai. Motivated by this, we consider next ui(k). 

Using the above notations, we see that 

Su (x,£,k) = S(x,£,T)(x))h(d(x,x ),x ,£o,k)ci(x)d 10 (x) G I°(SN;Yi), 

where C\(x) is a smooth non- vanishing function. Moreover, the operator 
Qk has the Schwartz kernel ( 1401 that away from the diagonal has the 
form 

Q k (x, f , x', £') = h(d(x, x'),x', k)8 + (x, £), 

where h is defined in (135ft . Thus, in (x, £, x', £') G Z \ L Q , the kernel of 
Qk has the form 

I ^ x ^ x, ^ fi) [h{d{x,x'),x',i',k)q{x,^e)]de modC°°(Z) 

where ip(x,£,,x',!;',9) is a non-degenerate phase function parameteriz- 
ing the LagrangianAp andg(x,£,#) G sn-i/2+(4n-2)/4-Jv/2^ R n xR n-i x 

M. N \ 0) has a non- vanishing principal symbol. 

Let us use in SN\rjo local coordinates S : {x, £) 1— > (sj(x, £)) •" 1 
having the property that if 7 X ^(R_) intersects the geodesic 70 then 
Si = si(a;, £) is the unique value such that 

7^(R_)n 7 o(K+) =7o(si), 
and s 2 (x,0 = d( l0 (si(x,O),x). By ]M Prop. 2.1], 

ut(k) = QkSuo(k) e r i --3(^;A 1 ,A 2 ) 
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and Ui(x, £, k) in (x,£) G SN \ r] Q has in the above local coordinates 
the form 



a(x, £, k) = h(s!(x, £), x , £ , /i(s 2 (x, 0, 7oOiO, 0)> C(z, 0> fc ) 

where £, #) is a non-generate phase function parametrizing the La- 
grangian manifold A 2 , ((x,£) = —jx,^(~ s 2(x,0) ^ s the direction of x 
from 7o(si) and p(x,£,,9) is a symbol with a non-vanishing principal 
symbol. Note that on A 2 \ Ai the principal symbol of a(x, £, k)p(x, £, 6) 
is non- vanishing on the conormal bundle of the submanifold 



where Si = Si(x, £), ^2 = s 2 (:r, £), C = CC^jOj an d &o(^,0 is non- 
vanishing and independent of k. 

Now we are ready prove unique solvability of the inverse problem. 

Proof of Theorem ll.2[ First we note that have found already the 
set 70 fl U . Thus we know the set W := SN \ (SM U r] ). By observing 
the singularities of u(k) at W, we can find the conormal bundle of 
the manifold K fl U. Thus by observing u(k) at W we can find all 
points (x,£) G W such that there is a broken geodesic from (x ,^ ) 
to (x,£) with a breaking point in M mt . Moreover, we can find the 
principal symbol of u(k) on N*K fl W in some local coordinates. By 
(|4~T1) . observing the asymptotics of the principal symbol on N*K fl W 
when k — > oo, we can find the function d(xo, 7o(si)) + <^(7o( s i) ; 
Si = Sx(x, £) on (sc, £) G W. Here 7o(si) G M m * is the point at which 
the broken geodesic from (xo,£o) to (x,£) breaks, that is, the broken 
geodesic changes its direction. 

Using the continuity of the geodesic flow, we can find all (x, £) G SN\ 
SM that are in the broken scattering relation with (x , £ ) and more- 
over, in such case we can find the broken geodesic distance d(x , 7o(«i))+ 
^(7o( s i)? This proves the result and even more: The singularities of 
the Schwartz kernel of the operator A determine the broken scattering 
relation R. □ 
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